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Abstract. Using numerical simulations of the magnetised solar photosphere and ra- 
diative diagnostics of the simulated photospheric models, we further analyse the phys- 
ical nature of magnetic photospheric intergranular vortices. We confirm the magnetic 
nature of the vortices and find that most MHD Umov-Poynting flux is produced by 
horizontal vortex motions in the magnetised intergranular lanes. In addition, we con- 
sider possible ways to directly observe photospheric magnetic vortices using spectropo- 
larimetry. Although horizontal plasma motions cannot be detected in the spectropolari- 
metric observations of solar disk centre, we find an observational signature of photo- 
spheric vortices in simulated observations of Stokes-V amplitude asymmetry close to 
the solar limb. Potential ways to find the vortices in the observations are discussed. 



1. Introduction 

The recent discovery of vortex motions in the solar photosphere has raised the question 
of the role of photospheric vortices in the energy transport from the solar interior to 
the outer atmosphere. Bonet et al. (2008) showed that Magnetic Bright Points (MBPs) 
are subject to vortex motions in the intergranular lanes. An analys is of small-scale inter- 
ranular vortex structures and their motions has been presented by'Varaas Dommguez et al.l_ 



(12011) . High-resolution multi-wavelength observations (Wedemeyer-Bohm & Rouppe van der Voort 



120091) suggest that the vor tices, correspond ing to MBP motions, coincide spatially 
with chromospheric swirls. IMoU et al] ( 20111) reported on the formation of small-scale 



intergranular v ortices in weak magnetic field simulations of the solai" photosphere. 
Shelvag et all (l2011al lbh have demonstrated that solar atmospheric vorticity genera- 



tion, which reveals itself as both MBP rotary motions and upper-atmospheric swirls, 
efficiently occurs only in intergranular magnetic field concentrations, and the vortic- 
ity generation in the upper solar photosphere is essentially a magnetic phenomenon. 
Vortex motions in small-scale magnetic field concentra t ions can generate a variety 
of mag neto-hydrodynamic waves (.Erdelyi & Fedunll2007l : IJess et al.ii2009, : fedun et al. 
201 Ibl) . which propagat e to the outer layers of solar atmosphere, where the energy can 



be dissipated (see e.g. Zaqarashvili et al.ll2011l) . This mechanism may be responsible 



for the transport of energy from the solar interior towards the higher layers of the solar 
atmosphere, including the corona. 

This paper is devoted to some further analysis of photospheric magnetic vortices. 
We demonstrate the generation of positive, directed outwards Umov-Poynting flux and 
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Figure 1. Horizontal slices of the modulus of the horizontal components of ve- 
locity (left panel) and vertical component of magnetic field (right panel) measured 
at the height y = 500 km above the average continuum formation level in the simu- 
lation domain. Regions with the vertical magnetic field strength larger than 400 G, 
bounded by contours on the left panel, are covered by small-scale vortex structures, 
visually different from granular outflows. 




also show that the vertical component of the Umov-Poynting flux vector is generated 
mainly by horizontal vortex motions. This finding is confirmed by the numerical simu- 
lation of an individual cylindrically symmetric magnetic flux tube embedded in a quiet 
convectively stable atmospheric model. We also perform detailed spectropolarimetric 
diagnostics with the Fel 630.2 nm line to explore possibilities for observing intergran- 
ular vortex motions close to the solar limb. 

The structure of this paper is organised as follows. In Section 2, we briefly describe 
our simulation setup and show the mechanisms to produce magnetic photospheric vor- 
ticity. We demonstrate the presence of Umov-Poynting flux connected to photospheric 
vortices in Section 3, and further support this evidence by performing a test flux tube 
simulation in Section 4, and discuss possibilities for observing photospheric vorticity 
in Section 5. Our conclusions are given in Section 6. 



2. Vorticity generation is the solar photosphere 



We use the MURaM code ( Vogler et al.ll2005l ) to undertake simulations of the magne- 
tised photosphere. The code is extensively used in the solar physics community and 



is des cribed in great detail in the literature (see e.g. iShelyag et al.ll2004l: iRernpel et al. 
20091) . Our simulation procedure is similar to that described in IShelvag et all (l2011ah . 



and here we provide only a brief description. The numerical domain is 12xl2x 1.4 Mm^ 
size, resolved by 480 x 480 x 100 grid cells. A uniform magnetic field of 200 G 
strength is introduced into a well-developed, non-magnetic photospheric convection 
snapshot, and a sequence of snapshots, containing the photospheric plasma parame- 
ters, is recorded. Here we use a well-developed snapshot of photospheric magneto- 
convection recorded after 1 hour of the physical simulation time, during which the 
magnetic field has collapsed into the intergranular lanes and the simulation has achieved 
statistical stability. 
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Figure 2. Dependence of the horizontally-averaged modulus of the vertical com- 
ponent of vorticity (left panel) and dependences of the Ti (dash-dot-dotted line), T2 
(dash-dotted line), (dashed line) and T4 (solid line) terms of vorticity equation 
(Eq. \T\ right panel) on height in the simulation domain. Zero height represents the 
approximate level of the photosphere. 

In Figure [T] the horizontal slices of the modulus of horizontal components of ve- 
locity V/, = ^v^ + vfj and the vertical component of the magnetic field By are shown in 
the left and right panels, respectively (note that in our notation x and z are the horizontal 
directions, while y is the vertical). The contours on the left panel bound the regions with 
B,, > 400G, corresponding to the magnetised intergranular lanes. It is evident from the 
figure that the intergranular vortex motions, visible as small-scale semi-circular struc- 
tures in the left panel of the plot, are bounded by the magnetic field contours, and, 
thus are clearly co-located with and abundantly cover the strong intergranular magnetic 
field regions. Also, there is clearly a difference in the geometry of the intergranular 
magnetic vortex flows and granular outflows. The granular flows have dark granular 
centres, corresponding to the upflows with small horizontal velocities, surrounded by 
stronger horizontal outflows, which do not show any small-scale structure. 

We analyse the small-scale vortex features in the magnetised intergranular lanes, 
using the vorticity equation, derived from the equations of magneto-hydrodynamics 
(MHD). The process of vorticity production was demonstrated by Shelyag et al. ( 201lj). 
Here we provide only a brief description. The MHD vorticity equation can be derived 
by taking the curl of the MHD momentum equation, separating the magnetic and non- 
magnetic terms and decomposing the magnetic term into two: 

J. ^2 ^3 ^4 

p^- -(ZrV)7-V-xV/?,-V-x[V/j^-(B-V)B] + ivx[(B-V)B]. (1) 

Dtp p p p 

Here, T\ is the vortex tilting term, T2 the baroclinic vorticity generation term, Tt, 
the magnetic baroclin ic term, and T/[ the magnetic field tension term. As shown by 
Shelyag et alj (1201 lal) . the latter term T4 is mainly responsible for vorticity production 
in the upper photosphere, while T2 generates vorticity deeper in the photosphere and 
in the convection zone. This is further demonstrated by the right panel of Figure |2j 
where the height dependences of the terms are plotted on the same scale. The figure 
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Figure 3. Horizontally-averaged dependences of the vertical component of 
Umov-Poynting flux vector (solid line) and the vertical component of the Umov- 
Poynting flux vector, generated by horizontal plasma motions (dash-dotted line), on 
depth. 



shows that T4 (solid line) becomes greater than T2 (dash-dotted line) at heights 0.2 Mm 
above the approximate continuum formation level, while T2 is about twice as large at 
the level of continuum formation 3^ = 0. The terms Ti and Tt, do not play a significant 
role in the vorticity production. It also should be noted that T2 and T4 act in an opposite 
way and suppress each other, as it is evident from the opposite signs of these terms in 
Equation [T] In the left panel of Figure |2j the dependence of the horizontally averaged 
modulus of the vertical component of vorticity is shown. It can be seen from the figure 
that, in accordance to the physical mechanism outlined above, the vorticity takes its 
local minimum aty = 0.2 Mm, takes its maximum, caused by the baroclinic term T2 at 
y = and increases towards the upper photosphere in the region y > 0.2 Mm. 



3. Umov-Poynting flux in the photospheric magnetic vortices 

The plasma motions in magnetic fields can generate Umov-Poynting flux. In magneto- 
hydrodynamics, the Umov-Poynting flux vector F is defined as 

F = B X (v X B) , (2) 

where v is the plasma velocity field, and B is the magnetic field. It is of interest to 
measure the vertical, (y) component of the Umov-Poynting flux vector (Fy), connected 
to the magnetic vortex motions in the solar photosphere, since it can be responsible 
for the energy transport from the solar interior to the outer solar atmosphere. Figure [3] 
shows the dependence of the horizontally-averaged vertical component of the Umov- 
Poynting flux vector on height in the simulated model. As the plot demonstrates, Fy is 
negative in the convection zone and in the low photosphere, while it becomes positive 
(directed outwards) in the upper part of the photosphere above y > 0.2 Mm, in the same 
range where the magnetic tension term is greater than the baroclinic hydrodynamic 
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Figure 4. Simulation of individual flux tube. Left image - magnetic field lines 
(red), velocity field lines (red-blue-green-violet) and total velocity at 1 Mm height 
(greyscale) in the domain. Right image - vertical component of Umov-Poynting flux. 



term, T4 > T2 (see Section 2). In our simulation we find the behaviour of the vertical 
component of the Umov-Poynting flux vector to be similar to the result demonstrated by 
ISteiner et al. (2 008) . However, it should be noted that Steiner et al. (2008 ) use a much 
weaker magnetic field in their simulation, and demonstrate a different mechanism of 
the Umov-Poynting flux generation. 

To confirm the connection between the photospheric magnetic vortices and the 
Umov-Poynting flux, we analyse the effect of horizontal motions on the vertical com- 
ponent of the Umov-Poynting flux by setting artificially the vertical component of ve- 
locity Vy to zero in Equation |2l The depth dependence of f v,v,=o computed for Vy = 
is shown in Figure [3] as the dash-dotted line. As can be seen, y^ -o is generally pos- 
itive, except from a small part of the domain deep in the convection zone, and f y,v,.=o 
is nearly equal to the full vertical component of the Umov-Poynting flux in the region 
y > 0.4 Mm. It is, therefore, clear that the positive vertical component of the Umov- 
Poynting flux in the upper photosphere is generated by horizontal vortex motions in the 
photospheric magnetic field concentrations. 



4. Vortex motions in individual magnetic flux tube embedded into a convectively 
stable solar atmospheric model 

To further establish a relation between photospheric vortex motions in the intergranular 
magnetic field concentrations and Umov-Poynting flux, we performed a test simula- 
tion of a cylindrically-symmetric magnetic flux tube embedded into a co nvectively sta 



ble m odel of the solar atmosphere, using the Sheffield Advanced Code (IShelyag et al. 
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20081) . In this simulation, the computational domain has a physical size of 1.2 Mm 
resolved by 100x100x196 grid points in the x, z and y directions respectively, with 
a lower boundary located at the ph otospheric level y - 0, and magnetic field con 



figuration as described in detail by iFedun et all ( 2011all3) . The expanding flux tube 



with the magnetic field strength of B = 1.2 kG at the footpoint (height y - 0.0 Mm) 
is in magneto-hydrostatic equilibrium with the non-magnetic ambient plasma outside. 
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Figure 5. Simulated Stokes-V amplitude asymmetry of the Fel 630.2 nm absorp- 
tion line at 60° angular distance to the north from the solar disk centre. 



We perturb the initially stable configuration with a velocity source located at height 
3^ = 0. 1 Mm in the photosphere and acting for the first 200 seconds of the simulation. 
The interaction of the source with the model creates a vortex with the horizontal veloc- 
ity amplitude v/, « 1 km s~^ at the footpoint of the magnetic f lux tube (see Figu relH this 
image was created with VAPOR 3D visualisation package ( Clyne et"aDl2007l ')). This 
vortex also generates the vertical Umov-Poynting flux, directed outwards to the upper 
solar atmosphere. The flux Fy is shown in the right panel of Figured Similarly to our 
previous analysis in Section 3, we compare the total vertical component of the Umov- 
Poynting flux with that generated by horizontal motions in the magnetic vortex. The 
two curves representing these are plotted as solid and dash-dotted lines, respectively. 
Since these are virtually indistinguishable, we conclude that the vertical component of 
Umov-Poynting flux is generated by horizontal vortex motions. It should be noted that 
only a qualitative comparison of this result with the Umov-Poynting flux amplitudes, 
provided in Section 3, is possible due to the different nature of these simulations. 



5. Radiative diagnostics of the photospheric magnetic vortices 

Intergranular vortices, being horizontal motions in the solar photosphere, are not easy 
to observe. The fast changing nature and small size of vortex structures, requires high 
spatial and temporal resolution. Such observations have been used to trac k spiral mo 



tions of photospheric features, including MBPs in the intergranular lanes (iBonet et al. 
llOOi). However, it is not possible to directly observe the velocity field of vortices in 
the upper photosphere using MBP tracking as the photospheric MBPs are localised in- 
tensity enhancements which are formed deep in the photosphere and do not fully cover 
the intergranular lanes. 

Therefore, we consider the possibility of directly observing photospheric magnetic 
vortex signatures using spectropolarimetry. Horizontal motions do not reveal them- 
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selves in spectropolarimetry at disk centre. However, the situation is different near the 
solar limb as the horizontal motions show a strong component along the line-of-sight. 
Also, due to the ma gnetic nature of the p hotospheric vorticity, the asymmetry of Stokes- 

V profiles (see e.g. IShelyag et al.ll2007l and the references therein), featuring in the re- 
gions with gradients of line-of-sight magnetic field and velocity due to layered structure 
of photospheric vortices (Shelyag et al. 2 011b ), is expected. We analyse this possibility 
by performing radiative transport diagnostics in the photospheric model, inclined by 
60 ° to the north from the solar disk centre, with the Fel 630.2 nm photospheric ab- 
sorption line. The synthetic Stokes- V amplitude asymmetry image, shown in Figure [H 
demonstrates strong signal in the intergranular lanes, as expected. Also, the Stokes- 

V asymmetry features in the image have s hapes closely resembling the cylinder-like, 
layered vortex structures, demonstrated by IShelyag et all fcOl lb ). One of the clearest 
features of this type is located at (3, 4) arcseconds in the image. The spatial resolution 
required to resolve this feature is estimated to be about 0. 1 arcsec in low light conditions 
due to limb darkening. Thus, future instruments for high-resolution spectropolarimet- 
ric solar observations with large apertures, such as the planned ATST, may be able to 
resolve these features in detail. 



6. Conclusions 

Using direct numerical modelling of solar photospheric plasma, we simulate the vortex 
motions in the magnetised intergranular lanes. Through two different approaches of 
simulations of magneto-convection, and direct modelling of the magnetic flux tube em- 
bedded into the solar atmosphere, we demonstrate that magnetised intergranular lanes 
generate significant amount of Umov-Poynting flux, directed outwards to the higher 
layers of the solar atmosphere, and show that the main part of this flux is produced by 
horizontal vortex motions. Finally, we demonstrate the possibility to directly observe 
photospheric vortices using high-resolution spectropolarimetric observations close to 
the solar limb. 
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